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Glass is a freezing phase of a deeply supercooled liquid. Despite its simple definition, the origin of glass 
forming ability (GFA) is still ambiguous, even for binary Cu-Zr alloys. Here, we directly study the stability of 
the supercooled Cu-Zr liquids where we find that Cu 64 Zr 36 at a supercooled temperature shows deeper 
undercoolability and longer persistence than other neighbouring compositions with an equivalent driving 
Gibbs free energy. This observation implies that the GFA of the Cu-Zr alloys is significantly affected by 
crystal-liquid interfacial free energy. In particular, the crystal-liquid interfacial free energy of Cu 64 Zr 36 in 
our measurement was higher than that of other neighbouring liquids and, coincidently a molecular 
dynamics simulation reveals a larger glass-glass interfacial energy value at this composition, which reflects 
more distinct configuration difference between liquid and crystal phase. The present results demonstrate 
that the higher crystal-liquid interfacial free energy is a prerequisite of good GFA of the Cu-Zr alloys. 
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Understanding glass-forming ability (GFA) 1 " 3 of metallic glasses (MGs) is a long-standing fundamental 
issue for basic science and their applications. With progress in strategies for alloying design, over 30 mm 
in diameter bulk metallic glasses (BMGs) prepared by directly copper mold cast technique were achieved 
in four alloy systems, i.e., 73 mm ZrCuAgAlBe 4 , 50 mm TiZrNiCuBe 5 , 40 mm PdNiCuP 6 , 35 mm 
LaCuAgAlCoNi 7 . However, the exact origin for the GFA in BMGs has not yet been completely understood 8-17 . 
Rather than the complexity in multi- component BMGs, binary Cu-Zr BMGs 11 (its critical diameter is larger than 
1 mm) have been intensively studied to figure out the GFA, because of its simplicity, and showing very sensitive 
dependence of composition on the critical size forming glassy phase, even narrow down to 1 at. %. 

Though enormous studies 11 " 17 , clear explanations have not been suggested for the origin of the good GFA for 
the Cu-Zr BMGs. For instance, in thermodynamic point of view, Kwon et al. 16 reported that the good GFA of 
Cu 64 Zr 36 was attributed to the relatively small difference of Gibbs free energy between amorphous and crystal 
phases. Wang et al. 11 pointed out that the GFA of Cu-Zr MGs was significantly influenced by the crystallization of 
competing phases. In addition, it was reported that the liquid structure of the Cu-Zr MGs could affect their 
thermal stability and the GFA; Yang et al. 15 addressed that a relatively high atomic-packing efficiency of short 
range order (SRO) on the Cu 64 Zr 36 BMG was attributed to the good GFA. On the other hand, in the kinetic point 
of view, Li et al. 12 demonstrated that the composition forming the BMGs on Cu-Zr alloys had higher density than 
other Cu-Zr compositions. This suggests denser metallic glasses have higher viscosity in the supercooled state, 
resulting in the good GFA. Moreover, Bendert et al. 17 suggested that the GFA of Cu-Zr alloys resulted from the 
different fragility at the supercooled liquid state and also near glass-transition temperature (T^), based on the 
measurement of thermal expansion coefficients. In other words, the stronger viscous behaviour of Cu-Zr liquid 
gives the better GFA. Nevertheless, the direct observation of such strong viscosity behaviour of the supercooled 
liquids forming Cu-Zr BMGs has not been reported yet. Therefore, the GFA is still ambiguous even for binary Cu- 
Zr MGs, although significant recent efforts have been devoted to study the GFA. 

One of determinant factors to facilitate the good GFA that has not been considered in detail is the stability of 
supercooled Cu-Zr liquids. Since glasses are a freezing phase of deeply supercooled liquid, the better stability 
should give the deeper undercooling, resulting in the vitrification of the supercooled liquid. Therefore, we should 
directly study the stability of the supercooled Cu-Zr liquids. We report crystal -liquid interfacial free energy and 
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glass-glass interfacial energy of Cu-Zr alloys, resulting in the stability 
of the supercooled Cu-Zr liquids using a containerless technique, i.e., 
electrostatic levitation (ESL), providing the deep supercooling, and 
using molecular dynamics (MD) simulation study. It is found that for 
a given composition the higher GFA, the higher undercoolability and 
the longer persistence of a supercooled liquid state. We provide 
compelling evidence to demonstrate that the interfacial free energy 
is a decisive factor to stabilize the supercooled liquid and thus affect 
the GFA, which is a novel approach to understand the GFA of the 
Cu-Zr alloys. 

Results 

Undercoolability and persistence of supercooled liquid. Figure la 
shows a representative heating and cooling curve of a Cu-Zr sample. 
The sample is heated over their liquidus temperature (T/) and cools 
down radiatively by turning off the heating lasers. Undercoolability is 
one of the signatures for the liquid stability, described by the degree 
of supercooling of the liquids, i.e., AT/T/ = (T/ — T r )/T/, where T/ and 
T r are the liquidus and recalescence temperatures, respectively. 
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Figure 1 | Time-Temperature curve of Cu 64 Zr 36 and undercoolability of 

Cuioo-xZr* (x = 35-54.3 at. %) alloys, (a) A typical Time-Temperature 
curve of the Cu 64 Zr 36 binary alloy by the electrostatic levitation (ESL) 
experiment, where T s is the solidus temperature, T/ is the liquidus 
temperature, and T r is the recalescence temperature, (b) Undercoolability, 
AT/T/ = (T/ — T r )/T/, of Qiioo-xZr* (x = 35-54.3 at. %) as a function of the 
Zr concentration. Three different liquidus temperatures of Cu-Zr alloys 
are used to determine the undercoolability, i.e., T/ from present ESL 
experiments (A), and from ref. 30 (•) and ref. 31 (■) phase diagrams. 
Meshed column is the result of the average liquidus temperature from the 
three liquidus T/. 



Figure lb shows the reduced undercooling values (AT/T/) of 
various Cu 10 o-xZr x (x = 35-54.3 at. %) alloy liquids. Interestingly, 
maxima are observed at Cu 64 Zr 36 (AT/T/ = 0.24 ± 0.01), Cu 56 Zr 44 
(AT/T/ = 0.30 ± 0.01) and Cu 5 oZr 50 (AT/T/ = 0.32 ± 0.01), 
regardless different liquidus temperatures. The maxima in 
undercooling are consistent with the critical thickness or the GFA 
on Cu-Zr MGs 12 , and thermal expansion coefficient maxima in Cu- 
Zr liquids 17 . The maxima in undercooling mean that the supercooled 
liquids of Cu 64 Zr 36 , Cu 56 Zr 44 and Cu 5 oZr 5 o have better stability, 
compared with neighbouring compositions. Based on the classical 
nucleation theory (CNT) 18 , the deeper undercooling reflects the 
higher nucleation barrier. Therefore, the liquids showing under- 
cooling maxima (i.e., better stability) should be expected to show 
longer persistence at a same degree of the supercooled temperature 
than the neighbouring composition liquids. 

To verify the stability of the supercooled liquid as a function of 
time, we carried out Time-Temperature-Transformation (TTT) 
experiments. In particular, we focus on the Cu 64 Zr 36 , since the origin 
of the higher GFA for the most studied composition is still unclear. 
Figure 2 shows the TTT-diagram of Qiioo-xZr* (x = 35-38.2 at. %) 
liquids. At the same reduced supercooling temperature, T/T/ = 
0.865, as an example, Cu 64 Zr 36 liquid indeed persists for much longer 
time 455 s in the supercooled liquid state than other neighbouring 
compositions, 124 s and 16 s for Cu 631 Zr 36 9 and Cu 6L8 Zr 38 2 , 
respectively. This experimental observation is consistent with the 
calculated TTT-diagram of Cu 64 Zr 36 (red) and Cu 618 Zr 38 2 (blue) 
from ref. 19, although the exact values are different. According to 
the CNT, the nucleation barrier (AG*) resisting the crystallization is 
determined by two factors, crystal-liquid interfacial free energy (er) 
and Gibbs free energy difference of liquid and crystal (AG V ), i.e., AG* 
= 167rcr 3 /3(AG v ) 2 18 . Here, the same reduced temperature is assumed 
to give the same driving force for the nucleation (in case of Turnbull's 
approximation, AG V = AHj(l — T/T/), where AHf is the fusion 
enthalpy). This assumption is reasonable because AHf should not 
be significantly different for the Cu-Zr liquids within such a small 
composition range (35-38.2 at. % of Zr). Accordingly, the longer 
persistence of the supercooled liquid at the same T/T/ indicates that 
the stability of the supercooled liquids is decisively determined by the 
interfacial free energy (a). 
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Figure 2 | Reduced Time-Temperature-Transformation (TTT) diagram 
of Cu 10 o-xZr x (x = 35-38.2 at. %) alloys. T/ at reduced temperature (T/T/) 
is the average value of the three different liquidus temperatures (T/: 
present, T z : ref. 30, T/: ref. 31). The blue and red lines are of calculated TTT 
curves from ref.20, and dots are of experimental results in present study 
using the ESL technique. The inset is magnified from the dots with dashed 
guide lines for eyes. 
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Interfacial free energy. Since the deeper undercooling and the 
longer persistence of the supercooled liquid Cu 64 Zr 36 imply higher 
interfacial free energy, to confirm it, we estimate the crystal-liquid 
interfacial free energy from the undercooling data using the CNT. To 
initiate crystallization reaction, at least one critical nucleus should be 
formed in the supercooled liquid. This requires that the number of 
the critical nucleus in the supercooled liquid should be greater than 
one with a given liquid volume V, a steady state nucleation rate I s , and 
time t at a certain temperature T 18 , i.e., 



I s (T)mVmt>l 



(1) 



Here, the steady state nucleation rate I s per unit volume at 
temperature T is given by 



I s = 



6 n* 2,3 k B TN A ( \b\i\ 



kt](T)1 3 \6nk B Tn* 



1/2 



exp 



k B T 



(2) 



where r\, X, n*, d/n, k B , N A and AG n * are the viscosity, the average 
atomic jump distance, the number of atoms in the critical nucleus, 
the Gibbs free energy difference between the initial and final phases 
per atom, the Boltzmann's constant, the Avogadro number, and the 
work (i.e., nucleation barrier) of critical cluster formation, 
respectively. The work AG„* is determined by the Gibbs free 
energy difference AG V = Sfi/v, the interfacial free energy (c), and 
the average atomic volume (v). As expected, it is found that the 
Cu 64 Zr 36 alloy, having the higher GFA 12 , exhibits the highest 
crystal -liquid interfacial free energy value in Figure 3, which is 
consistent with the undercoolability and TTT-diagram results in 
Figures lb and 2. It should be mentioned that the density and 
viscosity of Cu-Zr liquids measured in this study show weak 
composition dependence, which is reasonable with the narrow 
composition range. This is consistent with a previous report 17 that 
show no volume anomaly of Cu-Zr alloy liquids with compositions. 
Thus no composition effect of density and viscosity within the 
temperature range was found to estimate the interfacial free energy 
of Cu-Zr alloys. 

In addition, the glass-glass interfacial energies for Cu 10 o-xZr x (x = 
30-52 at. %) MGs were calculated at 800 K by molecular dynamics 
(MD) simulations in Figure 3b. It should be stressed that a prerequis- 
ite for calculating crystal -liquid interfacial free energies is the crystal 
structure, and the orientation of the crystal at the interface. 
Unfortunately, such information is not available as we change the 
composition of the alloys. And thus instead of directly calculating 
crystal-liquid interfacial free energies, we calculated the glass-glass 
interfacial energies of Cu-Zr MGs (see the details in Supporting 
Online Material). Since the difference of regularity or atomic packing 
density of short range order in bulk and surface of the glass is 
dependent of compositions, the interfacial energy between the glasses 
would show the dependence of composition. This gives the relation of 
short range order and GFA with composition. Although these data 
are not exactly the same as those for the crystal -liquid interfacial free 
energies, they still interestingly provide the same trends as the inter- 
facial free energies of crystal- supercooled liquid of Cu 10 o-xZr x (x = 
35-38.2 at. %) estimated from the undercooling data. 

Configurational differences causing interfacial free energy. Follow- 
ing the Turnbull's explanation of the stability of the supercooled 
liquid (or the high nucleation barrier) by crystal-liquid interface 
structure 20 , the higher interfacial free energy for Cu 64 Zr 36 alloy 
implies the larger difference of the configurational ordering 
between the supercooled liquid and the crystal. Unlike crystal- gas 
interface, atoms in the liquid near a crystal show ordering, called 
neg-entropic ordering. When the atomic configurational ordering 
at the interface is more different from that of the crystal, the larger 
interfacial free energy forms, and thus leads to the high nucleation 
barrier 21 ' 22 as follows. 
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Figure 3 | Interfacial free energies of Cu 10 o-xZr x (x = 35-38.2 at. %) 
alloys, (a) Estimated crystal-liquid interfacial free energies (cr) from the 
undercooling data in the ESL experiment using the classical nucleation 
theory (CNT) and Turnbull's coefficients (a) of Cui 0 o-xZr x (x = 35-38.2 
at. %) alloys, (b) Calculated glass-glass interfacial energy of Cu-Zr metallic 
glasses (MGs) at 800 K as a function of Zr composition together with the 
critical thickness for the glass formation in the Cu-Zr alloy system reported 
in ref. 12. 
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where N f is the number of atoms in the interface, N is the number of 
atoms in the crystal plane, AS con f ig (bulk) is the configurational 
entropy of the bulk crystal, AS con f ig (interface) is the configurational 
entropy of the interface, and ASf is the fusion entropy per atom. The 
higher Turnbull's coefficient (a) on Cu 64 Zr 36 (ot T = 0.519), as shown 
in Figure 3a, indicates that the configurational ordering between the 
crystal and the supercooled liquid of Cu 64 Zr 36 differs from that of 
neighbouring compositions. 

To further confirm the higher interfacial free energy for Cu 64 Zr 36 
alloy, we further carried out synchrotron X-ray diffraction mea- 
surements for Qiioo-xZr* (x = 35-38.2 at. %) MGs at ambient tem- 
perature. It is indeed found that the Cu 64 Zr 36 MG had the highest 
dense-packing of icosahedral short range ordering (ISRO) than other 
compositions, which is consistent with the recently reported results 
for Cu-Zr MGs 15 . Similar to other studies 23 , a high fraction of ico- 
sahedral-like clusters in Cu 64 Zr 36 MG was also detected after 
Voronoi tessellation analyses of the results obtained by reverse 
Monte Carlo (RMC) simulation to the experimental data 15 . Since 
the icosahedral short range order (ISRO) clearly differs from the 
short range orders in Cu-Zr crystal phases, the highest atomic 
packing on Cu 64 Zr 36 MG could result in the large difference of 
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configurational entropy between crystal and amorphous alloy. It is 
worth to note an opposite example which validates the conclusion; a 
liquid forming Ti-Zr-Ni icosahedral quasicrystal phase (/-phase) 
showed a lower interfacial free energy (er = 0.06 J/m 2 ) and ot T = 

0. 32 with small undercooling (AT/T/ = 0.1), since the both local 
orders of the liquid and /-phase are same local ordering, i.e, icosahed- 
ral order 24 ' 25 . 

The highest interfacial free energy on the Cu 64 Zr 36 liquid and glass 
links with the higher atomic-packing efficiency on Cu 64 Zr 36 BMG 15 , 

1. e., more difference of configurational ordering between the liquid 
and the crystal. Hence, the good GFA on Cu 64 Zr 36 comparing with 
neighbouring compositions is attributed to the stability of super- 
cooled liquid with the high interfacial free energy. Moreover, a recent 
simulation study 26 revealed the importance of crystal-liquid interface 
structure for the GFA; the poorer crystal-liquid interfacial structure, 
the larger difference of local ordering between liquid and crystal, 
leading to the better GFA for Cu 50 Zr 50 as compared with 
Ni 5 oAl 50 26 , which is consistent with this study. 

Discussion 

We carried out molecular dynamics (MD) simulation to calculate 
glass-glass interfacial energies for Cu 10 o-xZr x (x = 30-52 at. %) MGs 
at 800 K (Fig.3b). Interestingly, the glass-glass interfacial energies at 
800 K provide the same trends with the crystal- supercooled liquid 
interfacial free energy of Cu 10 o-xZr x (x = 35-38.2 at. %) estimated 
from undercooling data (Fig.3). Since the glass -glass interface is 
formed by welding two glasses, the interfacial energy reflects the 
different degree of randomness in short range order in bulk and 
interface of the glasses due to the more excessive volume on surface 
than in bulk of the glass. In other words, the surface has more defect- 
ive short range ordering 27 . Therefore, when the glass has the better 
regularity or the higher packing density of the ISRO in bulk, the 
larger difference of the excessive volume between the bulk and the 
surface is given, which causes the higher interfacial energy. 
Accordingly, the higher interfacial energy of Cu 64 Zr 36 compared 
with the neighbour compositions supports that the better regularity 
or higher packing density of ISRO results in the better GFA 15 . 
Moreover, other glass forming compositions, i.e., Cu 56 Zr 44 , and 
Cu 5 oZr 5 o also show higher glass-glass interfacial free energies. 

One would expect that the glass-glass interfacial energy would 
eventually diminish at high temperatures. One might also expect that 
such temperature is sensitive to alloy compositions, i.e., a higher 
interfacial energy usually corresponds to a higher diminishing tem- 
perature because the reduction of the interfacial energy at the high 
temperature is closely related to the glass relaxation time and the 
higher interfacial energy may be related to a higher viscosity. 
However, further simulation needs to be carried out to confirm this 
hypothesis. 

The C- curve in the TTT diagram is usually observed in bulk 
metallic glasses, due to the result of the competition of two character- 
istic times, i.e., x n and t d for nucleation and viscosity (or diffusion), 
respectively. Below a nose temperature in TTT diagram, t d for vis- 
cosity property governs the GFA. In contrast, x n for nucleation gov- 
erns the GFA above the nose temperature. Although electrostatic 
levitation can easily achieve deeper supercooling state than other 
methods, the cooling rate by radiation on Cu-Zr alloys is not fast 
enough to bypass the nose temperature in TTT diagram. This pre- 
vents from verifying the reported characteristic viscosity behaviour 
(i.e., fragility) 17 or the role of t d to explain the good GFA on Cu-Zr 
alloys. Instead, we can clearly see the characteristics of x n for nuc- 
leation with the higher undercoolability, the longer persistence of 
supercooled liquid state, and the higher interfacial free energy (i.e., 
nucleation barrier), which reflects the stability of the supercooled 
Cu-Zr liquids in figures 1, 2 and 3. Therefore, if a moderate cooling 
rate is given, the stability of supercooled Cu-Zr liquids should be 
kept, and the liquids can become glasses. 



In conclusion, we investigate the origin of the GFA on a Cu-Zr 
alloy system. We found that Cu 64 Zr 36 showed better undercoolability 
and stayed in supercooled liquid state for longer time, compared with 
neighbouring compositions. Under the same driving force, the 
longer persistence of the supercooled liquid on Cu 64 Zr 36 revealed 
that interfacial free energy was the key factor for the stability of the 
supercooled liquid. The highest interfacial free energy on Cu 64 Zr 36 
estimated from both undercooling experiment and simulation sup- 
ports this scenario, on comparing with the neighbouring composi- 
tions. Our findings could have an implication that the GFA for 
various compositions in a given alloy system is controlled by the 
interfacial free energy between crystal and supercooled alloy liquids. 
The higher interfacial free energy, the higher the GFA. This should 
certainly be validated in other MG systems. 

Methods 

Sample preparation. Spherical Cu 10 o-xZr x (x = 35-54.3 at. %) alloys were prepared 
by an arc-melting with high-purity copper (Cu, 99.99%) and zirconium (Zr, 99.995%) 
slugs under high-purity argon (Ar, 99.9999%) atmosphere. To avoid reaction with 
oxygen, a piece of Zr was used as an oxygen getter. The samples were melted at least 
four times for the homogeneity of each composition. 

Undercooling experiment and data acquisition. The levitated samples with 2- 
3 mm in diameter were melted under high vacuum condition in the range of 10~ 6 - 
10" 7 Torr by C0 2 heating lasers with triangle arrangement for symmetrical heating. 
Detailed description of the electrostatic levitation (ESL) equipment is given in refs.28 
and 29. The samples were heated over their liquidus temperature (T/) and cooled 
radiatively by turning off the heating lasers. The sample temperature was measured 
for every 6 ms by two infrared pyrometers with InGaAs detector operating at 
1.55 um wavelength (547-1973 K) and 1.6 um wavelength (773-2773 K). The 
spectral emissivity (E; ) of Cu-Zr alloy is determined with the known solidus 
temperature (T s ) by the phase diagram of the binary Cu-Zr alloys 30 . 

Synchrotron X-ray scattering experiment. For the structural study, the synchrotron 
X-ray scattering experiments were carried out in Hasylab in Germany 15 and in PAL in 
Republic of Korea. Related thermophysical parameters and MD simulations details 
are given in Supporting Online Material. 
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